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We report an ultrasonic study of the magneto-elastic coupling of the hydrogenated and deuterated 
(TMTTF)2PF6 organic salts. For both salts the temperature dependence of the longitudinal velocity 
along the c* axis displays a monotonic stiffening of the C33 compressibility modulus upon cooling. 
Below the characteristic temperature scale 40 K the modulus stiffening becomes markedly enhanced, 
in concomitance with the reduction of spin degrees of freedom previously seen in magnetic measure- 
ments as low dimensional precursors of the spin-Peierls transition. The magneto-elastic coupling 
appears to be much weaker in the hydrogenated salt due to the highly inhomogeneous elastic be- 
havior induced by the proximity of the charge ordering transition to the spin-Peierls phase. For the 
deuterated salt, an important anomaly in the ultrasound velocity is observed below the spin-Peierls 
transition temperature Tsp in agreement with scaling of the elastic deformation with the spin-Peierls 
order parameter. In spite of the weakly inhomogeneous character of the spin-Peierls phase transi- 
tion, the magnetic field dependence of Tsp is well captured with the mean-field prediction for the 
lattice distorted Heisenberg spin chain. 

PACS numbers: 75.30.Kz,75.40.Cx,75.80.-|-q 



I. INTRODUCTION 

The sulfur scries of quasi-onc-dimcnsional (ID) organic 
compounds (TMTTF)2X with the ccntro-synimetrical 
anions X ~ PFg, AsFg, SbFg... arc Mott insulators 
which develop a charge ordered (CO) state that 
is followed at lower temperature by either an antiferro- 
magnetic (AF) Neel or a lattice distorted spin-Peierls 
(SP) phase depending on the pressure applied, anion 
X substitution, and to some degree the deuteration of 
the methyl groups. In low pressure conditions, these 
(TMTTF)2A' salts stand out as model low-dimensional 
systems to study the competition between spin, charge, 
and lattice degrees of freedom. 

The (TMTTF)2PF6 compound is particularly well po- 
sitioned in the phase diagram allowing a systematic study 
of the interplay between differing degrees of freedom and 
the mechanism involved in the stabilization of a variety 
of phases. Indeed, this compound undergoes a SP transi- 
tion at ambientpressure followed by an AF ground state 
above 9 kbar [^-Q- As regards to spin degrees of free- 
dom under normal pressure conditions, the spin singlet 
state associated with the SP lattice distortion has been 
borne out by spin suscep tibility 0-0 > 

and NMR spin re- 
laxation rate |3, [13 [Hli whereas a sizable part of the 
magnetic field-temperature phase diagram for the PFg 
salt has been obtained by high-field NMR and magneti- 
zation studies In the hydrogenated (H12) PFk salt, 
the SP transition takes place at Tsp « 18 K [SJ?, Q 
and at 13 K in the deuterated (D12) compound [15|- In 



the II12-PF6 salt, the observation of X-ray 2fci? diffuse 
scattering indicates the presence below 60 K of one- 
dimensional lattice precursors of the SP transition [isj : 
these eventually lead to a reduction of spin degrees of 
freedom akin to a spin pseudogap [l^, as exhibited by 
EPR susceptibility 0,[i, and Q NMR data. 

If the role of lattice degrees of freedom in the stabiliza- 
tion of the CO ground state is relatively well documented 
in the PFg and AsFg salts [O, [3 , the situation is not 
so clear for the SP transition. Except for X-ray diffuse 
scattering for the lattice precursors below 60 K, which 
condense into Bragg spots at Tgp Il3t . and anomalies in 
the thermal expansivity data [171 llSj , the nature of the 
magneto-elastic coupling involved in the SP instability 
remains poorly understood in these compounds. As far 
as the charge degrees of freedom are concerned, they ap- 
pear largely decoupled from the spins in the non ordered 
temperature region. At Tsp however, a small dielectric 
anomaly has been observed in both II12 and D12 PFg 
salts, indicating a finite interaction between them as a 
result of three-dimensional ordering [l9j . 

In this paper, we use the ultrasound technique 
to investigate magneto-elastic coupling effects on the 
c* axis compressibility modulus of the II12 and D12 
(TMTTF)2PF6 salts. Our data reveal the highly inho- 
mogenous character of elastic properties of the H12 salt, 
compared to the D12 one, which is ascribed to the prox- 
imity of the CO and SP phases in the H12 salt. This inho- 
mogeneous phase yields in turn weaker magneto-elastic 
effects at low temperature when approaching the SP tran- 
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sition. Two anomalies are observed on the compressibil- 
ity modulus as the temperature is decreased: from 40 K, 
the modulus stiffening displays a marked enhancement, 
which is superimposed at T^p by a sharper increase, as a 
result of the SP order parameter. For the D12 salt, pre- 
cursor effects to the true 3D ordering appear at Tgp w 15 
K and are found to be inhomogeneous in character. The 
phase diagram obtained up to a 18 Tesla magnetic field is 
found to agree with previous measurements [l^ [l9j and 
the mean-field prediction for the lattice distorted Heisen- 
berg spin chain under field [20| . 



II. EXPERIMENT 

Hydrogenatcd (H12) (TMTTFjaPFg single crystals 
were grown from THF using the standard constant- 
current (low current density) electrochemical procedure. 
Starting from a dcuterated chlorokctone [l^, the syn- 
thesis of the fully dcuterated D12 was obtained us- 
ing usual procedures and the corresponding deuterium 
incorporation-levels is estimated by NMR to be 97.5 %. 
The crystals grow as needles oriented along the chain 
axis a and with two natural parallel faces perpendicular 
to the c* axis. We use a pulsed ultrasonic interferometer 
to measure the variation of the longitudinal acoustic ve- 
locity along c* relative to its value at To = 30 K, AVl /Vl 
= [VUT) - Vl(To)]/Vl(To). The acoustic pulses are gen- 
erated with LiNbOs piezoelectric transducers resonating 
at 30 MHz and odd overtones bonded to the crystals 
with silicone seal. Since the crystal structure is triclinic, 
the acoustic mode is characterized as quasi-longitudinal 
and is dominated by the C33 elastic constant or the c*- 
axis compressibility modulus. The ultrasonic technique 
is used in the transmission mode and, because the typical 
thickness of the crystal along the c*-axis is quite small 
around 0.4 mm, a CaF2 delay line must be used to sep- 
arate the first transmitted acoustic echo from the elec- 
tric pulse. Moreover, no transverse acoustic mode can 
be properly analyzed because of mode conversion and 
mode mixing at the different interfaces. The longitu- 
dinal mode can be measured because it has the largest 
velocity that permits its time separation from parasitic 
signals. Since the velocity is related to the mass density 
p and the C33 modulus through the relation C33 = pV^, 
the AVl/Vl data are directly the image of the relative 
variation of the compressibility modulus C33 if the den- 
sity variations can be ignored. Indeed, thermal expan- 
sion measurements predicts relative density varia- 
tion values one to two orders of magnitude smaller than 
the AVl/Vl data reported in this paper. Several crys- 
tals from different growth batches were used for these 
ultrasonic experiments. All crystals of the same com- 
pound showed similar overall temperature and magnetic 
field behaviors. The temperature was varied between 2 
and 140 K and a magnetic field up to 18 Tesla could be 
applied along the c*-axis. Unless otherwise stated in the 
discussion, all velocity data were obtained for a frequency 



around 100 MHz. 
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FIG. 1: (Color online) Temperature dependence of the rel- 
ative variation of the longitudinal velocity AVl/Vl- Upper 
panel PFe (H12): first (red) and third (blue) warming up to 
90 K; 24 hours relaxation at 3 K followed by heating up to 
140 K (green) and cooling to 2 K (black). Lower panel PFe 
(D12) and PFg (H12): cooling down from 140 K (blue and red) 
and warming to 140 K (green). The curves have been shifted 
vertically for comparison and the green curves are the same 
in both panels. Arrows indicate the transition temperatures 
defined in ref . [l7| . 



III. RESULTS AND DISCUSSION 

We present in Figure 1 the temperature dependence 
of AVl/Vl below 140 K for the two compounds. In the 
lower panel, the two top curves (blue and red) refer to 
the first cooling down of the D12 and the H12 crystals, re- 
spectively. At relatively high temperature the curves are 
noisy owing to the typical degradation of the interfaces 
between the organic crystal, the piezoelectric transducer 
and the CaF2 delay line upon cooling due to very different 
thermal expansion coefficients. Despite this degradation, 
we can identify two clear features for the D12 compound: 
beyond the typical velocity stiffening generally observed 
in solids as the temperature is decreased, there is an ex- 
tra enhancement in the stiffening below the temperature 
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scale Tint ~ 40 K noted in concordance with the charac- 
teristic feature seen in thermal expansion data [T7| . It is 
followed by another anomaly at much lower temperatures 
near Tsp. 

For the H12 compound, the same feature is present 
at the same Tj„(, but with a much smaller amplitude. 
However, it can be enhanced by warming the sample 
well above 40 K, that is in a temperature range where 
the thermal expansion coefficient along the c*-axis shows 
the strongest variation This process has yielded 

the green curve of the lower panel in Fig.l where the 
enhancement below Tint '~ 40 K sharpens significantly, 
contrasting with the virtual absence of anomaly at the 
true transition Tgp. 

Thus for the H12 compound, not only are we observing 
much smaller clastic anomalies, but it appears that the 
system enters a highly inhomogencous clastic phase be- 
low 100 K as it can be observed in the upper panel of Fig- 
ure 1 for different thermal cycles. When the temperature 
is maintained below 90 K, warming up the crystal several 
times produces a depression of AVl/Vl below Tint (red 
and blue curves, upper panel of Fig. 1). To retrieve the 
velocity enhancement below 40 K, we need first to relax 
the sample at 3 K during 24 hours, then to heat it up to 
140 K (green curve, upper panel of Fig. 1); subsequent 
cooling down to 2 K (black curve, upper panel of Fig. 1) 
reduces the enhancement below 40 K as observed during 
the first cooling down (red curve, lower panel). However, 
we observe then a separation of the two curves below 100 
K, which is indicative of an inhomogeneous state. Such a 
behavior is not found in the D12 compound over the same 
temperature range: repeated thermal cycles to 40 K yield 
always the same temperature dependence and, although 
only one other thermal cycling to 200 K (in a magnetic 
field of 14 Tesla) has been performed, no measurable in- 
homogeneity effects were detected. Below the crystalliza- 
tion transition near 200 K, the silicone bonds used with 
organic compounds arc highly stable and reproducible 
and, thus, the inhomogencous conditions observed here 
are intrinsic to the samples. They can be possibly linked 
with the closer proximity of the CO and SP orderings oc- 
curring at Tco = 63 K [ll,[2l| and Tsp = 16.5 K, respec- 
tively in the H12 compound, c omp ared to the Tco = 84 
K and Tsp = 13 K in the D12 [Hi. 

As to the nature of degrees of freedom that can be held 
responsible for this elastic behavior, those of the charge 
sector are unlikely to be involved, at the very least di- 
rectly since, according to electrical transport, charge ex- 
citations are known to be thermally activated up to rela- 
tively high temperature j2^ . However, a softening of the 
Young's modulus was reported around the CO transi- 
tion in other (TMTTF)2X salts, the amount of softening 
being 2% for AsFg, 0.5% for SbFg, and 3% for Re04 
(23j . This could be related to the charge disproportion- 
ation between adjacent TMTTF molecules to which the 
Young's modulus, dominated by the Cn elastic constant 
along the chain axis, could be more sensitive. For the 
c* axis, the C33 modulus rather probes the deformation 



between the organic and anion planes that remains es- 
sentially unaffected by the charge redistribution. This 
may explain why no similar elastic anomaly is detected 
on either D12 or H12 sample at the charge ordering scale, 
Tco, in the insulating state (see Figure 1). 

As for spin degrees of freedom, NMR and spin sus- 
ceptibility data display essentially gapless excitations in 
this sector 0, H, [l3| , down to the temperature scale of 60 
K where, according to early X-ray measurements, lattice 
precursors of the SP transition appear [Tsj . and below 
which spins and lattice become strongly coupled. Spins 
are known to couple easily to strain in various magnetic 
systems so as to induce a monotonic softening on the 
compressibility modulus whose temperature dependence 
correlates with the regular increase of magnetic suscep- 
tibility with temperature above the transition tempera- 
ture. For weakly coupled insulating spin chains for ex- 
ample (2^ . [2^ , the softening upon warming presents a 
maximum at a temperature comparable to the spin ex- 
change constant namely, T w J/ks, in agreement with 
the characteristic temperature for the maximum in the 
temperature dependent spin susceptibility [2^ . The situ- 
ation differs slightly for the (TMTTF)2PF6 salts consid- 
ered here. The insulating gap is actually relatively small 
[2^ [27I and the spins weakly localized [l|, so that the 
maximum in susceptibility is pushed above room tem- 
perature 0, 13 and so will be the maximum of the elastic 
softening. Thus the monotonic decay of spin susceptibil- 
ity appears consistent with the regular increase of sound 
velocity with decreasing temperature below 140 K (see 
Fig. 1). Moreover, the upturn in elastic stiffening below 
Tint ~ 40 K corresponds to a decline in magnetic suscep- 
tibility and NMR spin-lattice relaxation rate [4, j3, X, 9i] , 
which is indicative of the lost of spin degrees of free- 
dom due to lattice fluctuations [l^ [l^ . The interpreta- 
tion given here for Tint demarcates from the one given in 
the framework of thermal expansion experiments |l7l.[l8|. 
which associated Tint with different CO patterns or fer- 
roelectricity. 

Let us now examine the magneto-elastic coupling at 
lower temperature, namely in the vicinity of the SP tran- 
sition. Because of the highly inhomogeneous elastic char- 
acter of the H12 compound, we can only discuss the data 
for the D12 compound shown in Figure 2. Below 40 K, 
the AVl/Vl velocity data increases smoothly, presents a 
local maximum near Tgp w 15 K and then decreases ir- 
regularly down to a minimum value at Tgp = 13 K. In the 
SP ordered state, AVl /Vl starts increasing again with a 
temperature dependence reminiscent of the onset of an 
order parameter tied to a continuous phase transition. 

Although the velocity data do not indicate an inho- 
mogeneous character in the temperature range above 
Tgp , we have verified the stability of the elastic behavior 
against relaxation effects in the vicinity of the SP transi- 
tion. The curves shown in Figure 2 were all obtained in 
warming cycles from 2 K to 40 K, after different relax- 
ation periods at 3 K. Only small differences in amplitude 
arc found below ^ 20 K, but as regards the values of 
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source of the inhomogcneous SP phase. 
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FIG. 2: (Color online) Temperature dependence of the rel- 
ative variation of the longitudinal velocity AVl/Vl for the 
Di2 compound in the SP ground state below 30 K (warming 
cycle). The plots correspond to different relaxation periods 
at 3 K. Inset; data in zero field for the H12 compound where 
the small anomaly at Tsp is indicated by an arrow, the red 
curve being a guide to the eye. 




Tsp and Tgp both remain unaffected. The existence of 
hysteresis effects up to 20 K indicates that the system 
enters into some kind of inhomogcneous short-range or- 
dered phase, which is distinct from to true SP transition 
located at T^p = 13 K and where the uhrasound velocity 
exhibits its characteristic upturn. The low temperature 
limit attained by the velocity is found to be relaxation 
period dependent (see Fig. 2). This inhomogcneous char- 
acter linked to the SP phase was also observed in the mi- 
crowave dielectric measurements on the D12 salt but was 
absent for the H12 one [l^. In fact, the above elastic fea- 
tures tied to the SP ordering are much less pronounced 
for the H12 salt: a stiffening anomaly at Tsp w 16.5 K 
is scarcely visible through a change of slope in the ul- 
trasound velocity (see inset of Fig. 2). Because of the 
highly inhomogcneous elastic character below 40 K, this 
tiny anomaly is only observed during the first cooling 
down from 140 K and no hysteresis precursors appear to 
be observed above the transition at Tgp. It appears then 
that the elastic inhomogcneous character proceeds from 
two sources in these compounds: the outcome of a local 
electric polarization tied at high temperature to the spin 
degrees of freedom affecting the H12 salt over the full 
temperature range below 100 K due to the proximity of 
the CO and SP phases, and a short-range SP order below 
20 K preceding the true SP transition for the D12 com- 
pound. Neutron scattering experiments have revealed a 
weak three-dimensional (3D) SP amplitude in these salts 
(13 . ITsI likely originating from disorder in the 3D cou- 
pling between the chains; as this coupling proceeds from 
the PFg anion displacement, any inhomogeneities in the 
CO ordering will produce 3D disorder that may be the 
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FIG. 3: (Color online) Temperature dependence of the rel- 
ative variation of the longitudinal velocity AVl/Vl for the 
D12 compound in the SP ground state below 30 K at differ- 
ent magnetic field values. Arrows indicate the variation of 
temperature. The dotted line is a guide to the eye for the 
field variation of Tsp{H) and the dashed line indicates T*p 
below which an inhomogeneous SP phase appears. 

To go further into the analysis of the SP state, we 
present in Figure 3 the magnetic field dependence of the 
ultrasound velocity for the D12 compound below 40 K. 
Because of the above mentioned relaxation effects in this 
temperature range, all the curves have been obtained af- 
ter a stabilization period of 23 hours at 3 K before warm- 
ing the crystal up to 40 K (blue curve), followed by a 
cooling down to 2 K (red curve). The inhomogeneous 
character of the SP transition is clearly revealed upon 
warming (blue curves) and cooling (red curves) at all 
magnetic field values; the curves are clearly shifted from 
one another below 20 K and more particularly below Tgp, 
as indicated by the dashed line. When the magnetic field 
is progressively increased, Tsp decreases as expected for 
a SP transition [23|, but identically for both thermal cy- 
cles, which is consistent with the second order charac- 
ter of the SP transition. The situation is different for 
Tgp. which is field independent for the warming cycle, 
but decreases rapidly for the cooling one so as to merge 
rapidly with Tgp. Such an hysteresis loop around Tgp is 



5 



consistent with short range correlations or precursory ef- 
fects that are progressively accentuated by the magnetic 
field. Similar effects could not be found in the H12 com- 
pound for which the magneto-elastic coupling appears to 
be much smaller. 

As mentioned before, the AV^/^l data in the SP 
ground state mimics the temperature dependence of an 
order parameter. Indeed, elastic anomalies at a Peierls 
or spin-Peierls transition can be simply accounted for 
by a very simple Landau theory [28j . using two order 
parameters, the modulation (of wave vector q) propor- 
tional to the SP magnetic gap Ag and the uniform clas- 
tic deformation e. In addition to the SP condensate 
Fsp = a|A,|2 + b\Aq\^ and the elastic Fei = C§3e'^/2 en- 
ergies, the coupling energy Fc = h\Aq\'^e^ determines the 
type of elastic anomaly appearing at the phase transition, 
h being the coupling constant. Translational symmetry 
imposes the quadratic jA^p dependence; the term with 
r=2 is always allowed by symmetry, but r=l implies in- 
equivalence of |e| and -|e| elastic deformations. When 
minimizing the free energy with respect to the elastic de- 
formation e, one can show that the biquadratic {r=2) 
coupling term leads to a stiffening of the compressibility 
modulus C33 according to the expression, 

C33 = C3°3 + 2/l|A,|2, (1) 

where C33 is the bare modulus. 

On the other hand, to obtain a sudden softening at a 
transition (like at Tgp for instance) , a linear to quadratic 
coupling {r=l) is needed, which yields 

C33 = ^33 ^ (2) 

where b is the coefficient of the quartic term |Aq|^. 

We show in Figure 4 the square root of the velocity 
stiffening below Tgp, which is proportional to |Ag|, as 
a function of the reduced temperature T/Tsp. The am- 
plitude of the stiffening anomaly does seem to decrease 
monotonically with increasing field: according to Eq. 1, 
this means that either the SP gap |Ag| or the coupling 
constant h decreases with increasing field. For other 
SP materials, \Ag\ is not considered to be the singlet- 
triplet gap but a field independent energy required to 
break down a dimer in the SP ground state [2^, [2^ . If 
this is true for these organic salts, this means that the 
coupling constant is decreasing with field, h{H). Near 
the transition temperature the order parameter varies 
as Aq ~ (1 — T/Tsp)^, with a mean field exponent 
/3 ~ 0.49(2) in zero field, as indicated by the black 
curve in Figure 4. The exponent increases weakly with 
magnetic field. This mean-field value contrasts with the 
/? ~ 0.36 obtained for the dielectric constant enhance- 
ment Ae', measured on the hydrogenated compound, 
which was rather consistent with non classical exponent 
for a - one component - SP order parameter in three di- 
mensions [l9| . However, pre-transitional 3D fiuctuations 
are observed in the dielectric measurements, while they 
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FIG. 4: (Color online) Temperature dependence of the square 
root of the relative variation of the longitudinal velocity 
AVl/Vl as a function of the reduced temperature at fixed 
magnetic field values. The fuU fine is a fit to (1 - T/Tspf (/3 
= 0.49(2)). 

are practically absent in the clastic ones. Finally, the 
softening of the velocity data observed below Tgp down 
to Tsp is rather consistent with Eq. 2 and we must then 
consider the possibility that it may be dynamically en- 
hanced; as no frequency dependence effects have been 
investigated, we leave this as an open question. Since 
Tgp is not dependent on field, they can be considered as 
an energy precursor |A*| in an inhomogeneous SP phase. 




Tsp(H) /Tsp(O) 

FIG. 5: (Color online) Magnetic field dependence of the 
reduced SP transition temperature rsp(_H')/Tsp(0) for the 
deuterated (TMTTF)2PF6 salt: red triangles and blue 
squares are associated to different criteria defining Tsp{H); 
the continuous line is the best fit obtained from Ref.l9. 

The decay of the SP transition temperature under 
magnetic field is known to fit a quadratic variation in 
relatively low fields [l^, [25|, HO], which conforms with 
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the mean-field prediction for a quasi-one-dimensional SP 
system (2^, ■ From field dependent dielectric measure- 
ments at microwave frequencies [l9l |. it was established 
that the SP transition temperature for both the D12 and 
H12 salts follows the mean-field expression obtained by 
Cross m, 



Tsp(g) 
Tsp(O) 



1 - c 



U^rsp(o) 



(3) 



where the g-factor is fixed at 2 for (TMTTF)2PF6 
and the predicted value for the constant c is 14.4 [20|. Be- 
cause of the inhomogeneous character of the SP ground 
state shown in Figure 3, it is difficult to define only one 
criterion for Tsp{H). In Figure 5 we present the field 
dependence by using the 1^* (blue) and the 2"'' (red) 
temperature derivatives, the 2""^ derivative correspond- 
ing to the minimum of AVl/Vl (dotted line in Fig. 3). 
We have added the fit to Eq. 3 with c = 12.9 deduced 
from microwave dielectric data obtained on crystals of 
the same batch [l^- The agreement is good, but the dis- 
persion of data points yields a value c between 10 and 15 
due to the presence of the SP inhomogeneous phase. The 
fair agreement with the quasi-one-dimensional SP mean 
field prediction contrasts with the situation encountered 
in other SP systems like CuGcOa where pronounced de- 
viations from the predicted values are found [31 j . 



IV. CONCLUSION 

From ultrasonic measurements in hydrogenated and 
deuterated (TMTTF)2PF6 salts, we observed magneto- 



elastic coupling effects extending over a wide temperature 
range beyond the CO phase transition. In the hydro- 
genated salt, the relaxation dependent behavior of the 
compressibility modulus C33 is found to disclose highly 
inhomogeneous character of the non ordered SP phase, 
most likely due to the proximity of the SP and CO 
phase transitions in this material and to the amplitude of 
the spin-lattice coupling that is apparently much weaker 
compared to the deuterated one. For both compounds, 
there is an upturn in the modulus stiffening below 40 K, 
which is correlated to the loss of the spin degrees of free- 
dom observed in numerous magnetic experiments. An 
anomaly is clearly observed at the SP transition for the 
deuterated salt only. Although the SP phase appears 
only weakly inhomogeneous, the elastic anomaly at Tgp is 
consistent with a biquadratic coupling between the elas- 
tic deformation and the SP order parameter, a coupling 
that decreases in a magnetic field. Precursor softening 
effects at Tgp are rather explained by a linear-quadratic 
coupling. Despite the inhomogeneous phase, the mag- 
netic field dependence of the SP transition temperature 
agrees with the mean-field predictions. 
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